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ABSTRACT 


T t; in proponed Chat Jupiter's cloud bands represent 
Inri'o tumlo uonvoet Ion whom? chnrnclct 1 Is iletonulnod by the 
phase change of water at a level where the temperature is 
about 275K. It is argued that there are three important 
layers in the atmosphere: a tropopause layer where emission 
to space occurs, an intermediate layer between the tropopause 
and the water cloud base, and the deep layer below the 
water cloud. The barotroplc behavior observed in the visible 
clouds is seated in the tropopause layer where the stability 
is large. Horizontal potential temperature contrasts in the 
Intermediate layer are larger than vertical ones, are given 
by the difference between moist and dry adiabats, and pro- 
duce the observed thermal winds. Horizontal temperature 
contrasts in the deep layer are extremely small, and ordinary 
convection exists, with an unstable potential temperature 
gradient. The intermediate layer is stably stratified. 

It is argued that such a structure would spontaneously 
develop on a radiative time scale from an initial horizontally 
uniform convective state, and that the belt-zone spacing would 
be determined by the barotroplc stability criterion, as Inger- 
soll and Cuzzl (1969) suggested. All arguments are only semi- 
quantitative, since convection in a deep atmosphere or with 
phase changes is not understood. A major purpose is to point 
out that these ingredients are essential to Jovian meteorology. 



I. Introduction 

This paper outlines a theory of the Jovian general circula- 
tion. Latent heat release In water clouds and the lack of a lower 
boundary to the atmosphere are central to the theory. The 
discussion Is largely qualitative, because our understanding of 
the dynamical effects of phase changes and great depth is poor. 

A major purpose of the paper is to argue that further work In these 
areas Is necessary to understand Jupiter’s atmosphere, even to a 
qualitative first approximation. 

The arguments fall into three major categories, as follows: 

1) The observed wind speeds and belt spacings are consistent with 
an adiabatic structure varying between moist and dry in zones and 
belts. Ingersoll and Cuzzi (1969) showed that Peek's (1958) data 
on latitude variation of rotation period is consistent with winds 
satisfying the thermal wind equation (hydrostatic balance in the 
vertical and coriolis force balancing pressure gradient in the 
horizontal), with a fixed level of no motion at some depth and a 
fixed temperature difference between belts and zones. Barcilon 
and Qlerasch (1970) showed that the data is consistent with the 
fixed level being the water cloud base and with the fixed tempera- 
ture difference being the variation between the moist and dry 
adiabats, if the water abundance is given by the 

solar 0:H ratio. 

2) A hypothetical initial convective configuration, horizontally 
uniform on the belt-sone scale, would be unstable on a radiative 
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time scale to development of zones and belts. The reason is that 
the onset of a large circulation cell would quench moist convec- 
tion at downwelling locations, leading to radiative cooling of 
the layer between the level of emission to space (tropopause 
layer) and the water cloud base. Qierasch, Ingersoll, and 
Williams (1973) have shown that such a correlation between 
vertical displacement and radiative cooling can produce an in- 
stability leading to amplification of the initial displacement. 

o 

The time scale would be on the order of 10 s, consistent with 
observed variations in the large scale structure. 

3) The observed barotropic behavior of motions is consistent 
with a stable tropopause layer at the top of the atmosphere, 
and an Intermediate weakly stable layer between the water cloud 
base and the tropopause layer. It has been suspected for some 
time that the Jovian cloud motions are controlled by two- 
dimensional vorticity conservation (Ingersoll and Cuzzl, 1 969). 
Striking qualitative confirmation now exists (Williams, 1975 i 
Maxworthy, 1975 > Ingersoll, 1973). We argue that such behavior 
is to be expected if the Rossby number e satisfies e << 1, if in 

p 

the tropopause layer e Ri >> 1 (where Ri is the Richardson number), 

p 

and if e Ri << 1 < Ri in the intermediate layer. We also argue that 
the moisture-related Instability should lead to this structure. 

These points are not treated in this order, or even systema- 
tically one at a time. It is more convenient, following the 
equations of motion, to treat short time scales first (Section II), 
and then long time scales (Section III). In Section IV the 
uniqueness of Jupiter and applications to the other outer planets 
are discussed. 
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An excellent and comprehensive review paper has just been 
written by Stone (1975) and It is not necessary for us to discuss 
previous theoretical modeling here in any detail. Realistic treatment 
of vortical structure has not been attempted. An artificial lower 
boundary has generally been assumed, removing the coupling to deep 
layers which exists on the real planet. The consequences are 
both thermal and dynamical. Deep layers are alnost certainly 
very close to horizontally isothermal (Bogart and Qierasch, 1976) 
and a boundary removes the constraint that should be imposed on 
the atmosphere. A boundary has the dynamical effect of permitting 
horizontal pressure gradients to exist at the base. We comment 
further on this in Section II. The other deficiency in previous 
work has been neglect of latent heat release. On Jupiter this is 
one of the largest potential sources for horizontal temperature 
variation. 

An exception is the work by Barcilon and Glerasch 
(1970). An attempt was made to treat both latent heat and 
coupling to deep levels. However, the discussion was based on 
a highly Idealized axisymmetrlc steady flow, vith a frictionally 
controlled meridional circulation. No attempt was made to explain 
self-con3istentiy the origin of the flow or its scale. 

We shall make several basic assumptions throughout this 
paper. One is that the abundance of water on Jupiter is given 
to correct order of magnitude by the solar 0:H ratio. For 
levels of cloud condensation on different planets we refer to 
Weidenschllling and Lewis (1973). We assume that condensation 
of water leadB to precipitation, so that there is differentiation 
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between dry and moist gas. This Is discussed by Rossow and 
Qierasch (1976). We assume the ideal gas equation of state. We 
use locally cartesian geometry (a B-plane), since the scale of 
the belts and zones is much less than the radiUB of the planet. 
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II. Short Time Scales 

In this section we discuss motions whose time scales are 
on the order of L/U, where U is the order of the flow speed 
and h is the characteristic horizontal scale. The important 
motions in this category are the shear instabilities at cloud 
level, and other instabilities that might exist at other levels. 

The length scale L * D/2tt, where D is the spacing between 
zones. With D » 10^ cm and U ■ 25 ms"* 1 (consistent with a dif- 
ferential rotation period of about one minute) we obtain the 
time scale L/U ^ 10^s. We assume that friction, radiative 
heating, and eddy diffusion are negligible on this time scale. 

We shall also assume that the depth of the flow is much less 
than L, so that the hydrostatic approximation is valid. The 
equations in 'pressure coordinates are 


Du 

Dt 


fv 



0 


( 1 ) 


Dv 

Dt 


+ 


fu 


. 3<J> 
+ 37 


( 2 ) 


3$ 

3p = 




( 5 ) 


a 
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a v , aw 
57 + 3p 


0, 


C*D 


PJS> „ 

Dt 


o, 


(5) 


where 
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. Jfci 
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T f k - Jrr » f * f rt + 3y. 


and the notation Is 


x, y 
u, v 
P» w 
t 

D/Dt 

* 

Y, R 
T, © 

Po 

V 0 


■ eastward, northward directions, 

* x, y velocities, 

■ pressure, Dp/Dt, 

■ time, 

« 3/3t + u3/3x + v 3/3y + w 3/3p, 

» geopotential height, 

» specific heat ratio, gas constant, 

« temperature, potential temperature, 

* reference pressure (water cloud base) 

■ coriolls parameter, and its gradient. 


Let © and 4> be written as a static part plus a part associated 
with dynamics: 

(6) 


© = 6 s (p) + 6(x,y,p,t) 


* = 4> s (p) + 4>(x,y,p,t) 


Let the amplitude A6 be associated with 0. Define 


e = 


U 

f o L * 


« _ £o d0S n 

A0 dp * f 


_ 6L 


Note that S(p) is a function of depth; e and B are not. 
Nondimensionalize as follows: 

© by A0, x, y, p by L, L,p Q 

<J> by RA0, t by U/L, 

u, v by RA0/f Q L = U, u by p Q U/L . 


(7) 


( 8 ) 


( 9 ) 
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Equations (1) - (5) become, 

* ' (1 * Bj,)v + 5x 

c BF + (l + Bv )u + ff 


in terms of dimensionless 

variables 

- 0 

(10) 

- 0 

(11) 

SC 

° HI 
C 

1 

■ 

(12) 


du 

3x 


+ 



dU 

5p 


!&♦* 


■ 0 


« 0 


(13) 

(14) 


We proceed by expanding in e, as is usual in developing 

_n 

quasi-geostrophic theory. Observationally , U •v* 25 ms , 

0 _ h — 1 

L ^ 2 x 10° cm, and f * 2 x 10” s . , so that e^i/16. Let 
b « B/r ; observationally, b is order unity. Write all dependent 
variables as a power series In e; for example, v * v^°^ + ev^^ 

+ ... . To leading order we obtain from (10) - (13) 


„<o> M <0> . „ 
v + 3x °» 

(15) 

u (o) + i± ^ m 0 

U + 3y U > 

(16) 

34» 6 

8p * - pl-K ’ 

(17) 

(o) 


_ o 

3p °* 

(18) 


We do not know how to write (14) to lowest order, since we do 
not know the order of S. 
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Going to first order In the dynamical equations (10) and 
(11), using (13), (15), and (16), we obtain the equation for 
the vertical vorticity 


0^°) (a/ 0 * 

Dt \3x 





(19) 


where D^°VDt is evaluated with u^°\ v^ 0 ^. We have assumed 
the upper boundary condition *► 0 as p *► Oj together 

with (18) this implies that u) (o) » 0. We justify this 
boundary condition below. Finally, (l 1 !) gives to leading 
order 


D (o) 

"BtT - 



+ 


eSw 


( 1 ) 


0 . 


( 20 ) 


This completes the development of the equations, and we proceed 
to discuss application to Jupiter. 

Rather than beginr' g by attempting to estimate the un- 
certain quantity S, let us first list the three qualitatively 
different possibilities for the system of equations (19) -ini (20). 


1. eS « 0(1). Equations (19) and (20) are coupled. This is 
the standard quasi-geostrophic approximation which successfully 
predicts the major features of mid-latitude circulation on Earth. 
The thermal field influences the vorticity by affecting the 
height of columns of fluid between isentropic surfaces, and 
vorticity advectlon affects the thermal field by causing vertical 
motion. When horizontal temperature gradients are present, 
baroclinic instabilities can develop, and the length scale of 
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the Instabilities Is Internally determined to be precisely 
that which couples (19) and (20): 


r „ . Po nd V dp 


“ 0 ( 1 ). 


The L determined this way Is called the rtossby radius of 
deformation. 

2. eS >> 1. Equation (20) predicts » 0, and (19) reduces 

to vorticlty conservation. The thermal field Is secondary. 

♦ 

3. eS << 1. Equation (20) predicts D (o V°VDt - 0, (19) 

I 

leads to vertical motion through vorticlty advection, but the 
vorticlty and the vertical motion field are secondary. 


We propose that there Is an abrupt transition In the 
Jovian atmosphere between eS >> 1 at tropopause level and 
eS « 1 at deeper levels, and that the complicated intermediate 
regime typical of the Earth's atmosphere does not occur on Jupiter. 
We believe the observed barotropic behavior of the cloud motions 
supports this proposition. Furthermore, if one accepts the 
water cloud hypothesis for the origin of the horizontal tempera- 
ture contrasts as, it can be argued quite reasonably that the 

transition will exist. Baroilon and Gierasch (1970)show that 

_2 

latent heat can lead to A0/O ^ 10 . The atmosphere divides 

into three regions: 

1. Tropopause layer. This is a layer about one scale height 
deep where emission to space takes place. Radiation stabilizes 
this layer, so that 


Then eS can be written either of two ways: 


_ p o 0 s 

- eS ' c — So 



2 

jnfr * NT 


» 


where we used p Q /p 14, R « 4 x 10^ erg K” 1 gm*’ 1 , 0 g ■ 250 K, 
f - 1.7 x 10“ J *5"' l » and L - 2 x 10 8 cm.“ 

2. Intermediate layer. This is the region between the tropo- 
pause layer and the base of the water cloud. Since the maximum 
potential temperature contrast available here is A0, we expect 
dS , 'dp £ - A0/p o . Then 


-eS 


* E 


3. Deep layer. This is below the water cloud. Here there 
are no latent heat sources, and radiation is not effective. 

We expect ordinary convection to exist. Mixing length theory, 
for example, predicts 1 >> S > 0. 

This completes the general outline of the structure we 
propose. Figure 1 schematically illustrates the layers. We 
proceed now to a few specific remarks concerning the behavior 
of solutions. Discussion of the 3low processes which are 


* All quantities in this paragraph are dimensional, since we 
are evaluating e and S, not using the equations of motion. For 
a summary of numerical values and reference to their sources, 
see Table 3 in Section IV. 
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responsible for the genesis and maintenance of the banded 
structure is left for Section III. Here we deal only with 
questions related to short term stability. 

1. Conditions for a steady flow. The geostrophic and hydro- 
static relations (15)-(17) can be used to solve for and 

v (o) in terms of If 0^ is to be steady in the inter- 

mediate layer, we must have 

(o) 30^°^ , (o) 30^°^ 

u' ' + v x/ ■ o . 

and therefore 

3 * 9 ^ 39 ^ 3 2 9 ^ o) 39 ^ 

9y3p 3x 3x3p 3y 0 * 

From this it follows that the quotient (30^°V3y)/(30^° Vsx) 
is a function of x and y alone, and therefore that 

9 (o) » F(p)G(x,y ) , (21) 

that is, the isentropes are similar in shape at different 
heights . 

2. Stability of the intermediate layer. It is rather difficult 

to imagine how the structure (21) can arise and then be maintained. 
We propose here that it obtains automatically because F(p) is 
nearly constant; that is, horizontal contrasts in the potential 
temperature are larger than vertical ones in tne intermediate 
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layer. The reasoning is as follows. In Section III we present 
a mechanism for generating horizontal temperature gradients on 
the belt-zone scale. Vertical gradients In the potential tempera- 
ture do not arise directly. But it is well-known that in the 
presence of a horizontal gradient, instabilities can set in with 
the effect of converting some of the potential energy available 
because of the horizontal contrast into stable layering in the 
vertical. In particular, one class of instability is Inertial 
instability, which Stone (1966, 1967, 1971) has discussed in 
connection with Jupiter. These are small scale instabilities 
which can exist locally, Independent of the nature of boundary 
conditions (McIntyre, 1970). 

Inertial Instabilities grow rapidly (with a time scale 

—1 

on the order of f ) when the Richardson number, Ri, is less 
than unity. Ri is given, approximately, by 

d0 s ! 

M-vSit * § ; 

u 2 e 

We propose that these instabilities operate to adjust the 
stability so that Ri > 1, and as a result, so that S ^ e. It is 
then clear from the definition of S that the ratio of vertical to 

I 

horizontal potential temperature contrasts is small. , 

I 

The real question is of course whether other processes or 
other instabilities act to stabilize the intermediate layer still 
further. No other instabilities are known at present to exist i 

under these circumstances (e << 1, Ri > 1) in an unbounded atmos- <• 

phere, but this is clearly an Important question for future work. 
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3. Stability of the upper layer. Plow in t Is layer is 
governed by the barotroplc vortlcity equation, and the stability 
criteria should be the familiar ones for this kind of flow. 

For example, for a purely zonal mean flow the criterion is 


b 


9 2 u (o) 

— jp- 


> 0 , 


( 22 ) 


and this is one of the results supporting the proposition we 
have made. 

Vertical motion and interaction with deep layers. In the 
event that barotroplc instability sets in at the tropopause 
layer, the resulting mixing could lead to unbalanced pressure 
gradients which drive horizontal convergence. The horizontal 
"friction" caused this way would represent a breakdown in our 
initial assumptions. 


The resulting adjustments would change the thermal structure 
in the intermediate layer to make it consistent with the new 
flow pattern at tropopause level. Thid process is extremely 
complicated, and a detailed numerical calculation would probably 
be required to investigate it. One question of interest, how- 
ever, is the influence of such adjustments on the deep layer. 

Here we shall show qualitatively that the penetration depth 
is very large. 

Combining (19) and (20) gives 


,(o> 

"dF" 



3 

3p 


1-K 



* 


f 


"Sx ! 


+ 


+ by - 


S 


(23) 


in the deep layer is very small, ao that (23) 


Aaaume that 9' ' 

can be linearized * Let it also be proportional to exp(lot + iXy 
+ Ikx ) , i.e., Fourier analyze the disturbance. Then (23) gives 



Since we do not know the sign of the factor multiplying 9 in 
(2 1 !), we do not know if the solution decays with depth or not. 

In either case, the qualitative nature of the solution can be 

*1 wm • 

obtained by assuming that S/p is constant. The solution 
then gives 

■ <0 [Hfa •>' • ¥ >1' »->.’] • <»> 


+ X 2 - 


where A is an arbitrary constant, and the choice of sign and 
function depends on the sign of the factor on 9 in (2*0. We 
have assumed that the correct boundary condition is = 0 
at some deep level p = p^. We see that as long as S << 1 in 
the deep layer, the important results are independent of all 
other parameters: 
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a) The penetration depth for la large; If the stability 

S Is small enough, the depth Is determined by whatever interface 

I 

or phase change exists to prevent motion. 

b) To leading order, the boundary condition Imposed on the 

intermediate layer by the deep layer is * 0. The vertical 

motion is unconstrained. 
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III. Long Time Scales 

Here we discuss slow processes whose time scale Is governed 
by the thermal inertia of the atmosphere. This time scale Is 

c P-AO fl 

fc * * 1.3 x l(Ts 4.1 yr.), 

O 1 1 Q 

where we took c p « 1,3 x 10 erg gm" A K“ , p Q * 7 x 10 erg cm*’' 3 , 

A0 * 1.7 K» g « 2400 cm s~ 2 , and P ■ 5 x lO^erg cm“ 2 s” 1 . This Is 
the time required for the energy flux characterizing the system 
to lead to temperature changes on the order of A0. The Important 
slow processes are, of course, the genesis and maintenance of 
the belts, zones and other permanent features of the general 
circulation. 

Consider first a hypothetical initial state without belts 
or zones. The internal heat Is carried by small scale convection. 
At very deep levels, horizontal temperature gradients are con- 
strained to he extremely small because the superadiabaticity is 
small. Bogart and Qierasch (1976) have shown that for ordinary 
mixing length theory with rotation neglected, horizontal gradients 
on the order of the superadiabaticity lead to efficient horizontal 
redistribution of heat, so that even when an uneven heat flux is 
applied to the top or bottom of a deep convecting layer, horizontal 
temperature gradients remain small. It is unlikely that the 
influence of rotation in Jupiter is strong enough to change this 
conclusion. 
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The uneven heat flux applied at the top of the Jovian atmos- 
phere by the latitude variation of the Incident solar radiation 
therefore leads to heat redistribution at great depth, and the 
extremely small horizontal temperature contrasts necessary to 
accomplish the redistribution explain the observed lack of latitude 
variation in emission to space (Ingersoll et al., 1975b). 

An important point is that the internal heat is large enough 
so that the convective heat flux exists at all latitudes. It is 
minimum at the equator, where the largest fraction of the emission 
to space is supplied by solar heat. But the ratio of total 
emission to total solar energy absorption would need to be less 
than about 1.3 for convection to cease entirely at the equator, 
as opposed to the observed 1.9 (Ingersoll et al., 1975a). 

A few details of this hypothetical convective initial state 

are: 

1) The moist and dry adiabatic lapse rates differ 
appreciably only In a thin layer near the water cloud (Barcilon 
and Qierasch, 1970). However, the latent heat may influence 
the nature of convection drastically by introducing larger 
temperature perturbations, as Is the case on Earth. 

2) Radiation is important only in the tropopause layer, 
where cooling to space balances the convective flux from below. 

3) Neglecting rotation and latent heat, and assuming 
that the mixing length is one pressure scale height, properties 
of the convection are: 


t 
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RT 
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-1/y 


(27) 

( 28 ) 

(29) 

(30) 

(3D 


where F la the convective heat flux, w is the convective velocity, 
p la density, a subscript zero denotes reference level (water 
cloud base), H is scale height, g is gravity, K is eddy diffusivity, 
and T' is temperature perturbation. Numerical values will be 
of interest, and a few are given in Table 1. 


Table 1 

Mixing-length convection properties 

These are evaluated from (?7)-(31), with R = 3.7 x 10^ erg gnT^K” 1 , 
T “ 275 K, g B 2*100 cm s” 2 , ,p 0 = 7 x 10 6 erg cm~3, y ■ i.i| t and 

F = 5 x 10° erg cm" 2 s"l. 


p/p o w(cm s H(cm) K(cm 2 s" 1 ) t(s) T'/T 


0.1 

4.7 

X 

10? 

2.0 

X 

log 

9.4 

X 

iog 

4.2 

X 

10^ 

4.1 

X 

10 

0.3 

3.6 

X 

1°2 

2.8 

X 

10 6 

1.0 

X 

10 Q 

7.8 

X 

iog 

1.8 

X 

10 

1.0 

2.7 

X 

io; 

3.9 

X 

iog 

1.1 

X 

10? 

1.4 

X 

10? 

7.2 

X 

10 

3.0 

2.1 

X 

io; 

5.3 

X 

10 6 

1.1 

X 

10 Q 

2.5 

X 

10? 

3.1 

X 

10 

10.0 

1.6 

X 

10* 

7.5 

X 

10° 

1.2 

X 

10 9 

4.7 

X 

10^ 

1.2 

X 

10 


Notice that the time scale t is short enough so that rotation is 
probably not a dominant effect except at levels deeper than tfiese. 
The values of T'/T show that the superadlabaticity is extremely 
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small. A time scale for horizontal mixing over a distance 
L ■ 2 x 10® cm would be t^ L 2 /K. Evaluated at p/p Q - 1 , we 
find tj ^ 3.8 x 10^a l year; notice that this Is slightly 

shorter than the thermal time scale evaluated above. 

4) In reality, latent heat effects would lead to a mean 

thermal structure somewhere between the moist and dry adlabat3, 

If observations in the terrestrial atmosphere provide an accurate 
guide. The reason Is that a moist adiabatic structure is actually 
stable because rising bubbles of moist material always entrain a 
certain amount of dry gas during their ascent through dry layers . 
For discussion, see e.g., Ooyama (1971). Unfortunately, our 

t 

understanding of these processes is not complete enough to extend 
ideas to Jupiter and make predictions. Furthermore, molecular 
weight differences are much greater on Jupiter than on Earth, 
and are more important. 

This completes our discussion of the hypothetical convective 
basic state. The important feature of the configuration is that 
it has adjusted itself to be only marginally unstable against 
perturbations on a dynamical (convective) time scale. We now pro- 
ceed to argue that it would be unstable on a radiative time scale 
against certain large scale flows. We shall assume, although it 
is not crucial, that the basic state thermal structure is approxi- 
mately given by the moist adiabat. 

Consider the sequence of events: 

1) A weak large scale flow develops, with downward motion at 


certain locations. 
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2) At these locations, dry gas now exists at the level where 
moist convection was previously occurring. 

3) Convection ceases at these locations because the mean thermal 
structure is uubadiubutle with respect to dry convection. 

*0 Radiative cooling continues at the tropopause level and there- 
fore in the absence of convective heating from below, the locations 
with downwelling begin to cool, at all levels between the tropo- 
pause and the water cloud. 

5) At locations with upwelling, the mean thermal structure re- 
mains unchanged. It is already close to the moist adiabat and 
cannot be heated further by convection. 

On the average, then, there is a positive correlation between 
radiative heating rate and vertical displacement. Qierasch, Inger- 
soll, and Williams (1973) have shown that this can lead to in- 
stability, Gierasch (1973) shows that in equatorial regions, 
zonal symmetry is favored. For details, the reader is referred to 
these papers. We now list a few comments on the proposed instability. 
1) The Gierasch, Ingersoll and Williams discussion assumed that 
vertical motion leads to increased infrared blanketing, and there- 
fore radiative heating. An assumed static stability (sub- 
adiabatic structure) is necessary. The formulation is therefore 
self-consistent only if the circulation is seated at high levels 
in the atmosphere so that the motionless basic state is sub- 
adiabatic. Here we are proposing that the circulation extends 
down to the water cloud level, and we achieve a subadiabatic 
structure on the average because the moist adiabat is stable 
against dry convection. 



21 



2) The amplitude of the perturbation would be limited by the 
fact that AO could not exceed the difference betwet the moist 
and dry adlnbats. For Jupiter* assuming a solar abundance of 
oxygon, this Is AO *** 1.7 K (m?o {tout Ion IV, Table 2 for references 
and discussion). If we form a thermal wind equation by eliminating 
* between (2) and (3), and then integrate from p Q to p^, where 
P 1 is the observed cloud layer, we obtain 


u(p i) 


b agi 1 ~(Pi/p 0 ) 

f 0 3y k 


(32) 


where we assumed u(p Q ) « 0. Numerical values are presented in 
Section IV; for now we wish to point out how the scale of latitude 
variations might be determined. The radiative instability 
analyses showed that the smallest scales are favored. However, 
small scale belts and zones would be unstable on a short time 
scale to barotroplc instability. The criterion is 


3 - 



* 0 . 


(33) 

\ 


If we assume that A0 does indeed always reach, but never exceed, 
the moist-dry adiabat difference, and if we assume that uCp^y) 
is a smoothly varying function of y (as is suggested by the in- 
stability analyses in the linear regime), then we see that a 
scale is. uniquely determined between (32) and (33). This result 
is consistent with Ingersoll and Cuzzl's (1969) discussion of 
Peek's (1958) data. The reasoning is obviously far from complete 
without detailed solutions, but nevertheless is suggestive. 


I 
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3) The amplitude of the vertical motion field would be limited 
by the available heat flux. The heat flux Is given by P pc wT', 

r 

where w Is the vertical velocity magnitude and T’ Is the tempera- 
ture difference between rising and sinking gas. With 

T. .(£|‘a6 . 

and AO ■ 1.7 K, we reach the numerical values listed in Table 2. 

We have used (29) and (30) to estimate diffusion coefficients arid 
characteristic times for each level, although the Interpretation 
is clearly different from the mixing length case, since a single 
large cell In the vertical Is assumed. The Important point is 
that motions are much slower in the present case, because a much 
larger temperature difference is introduced by the latent heat. 

Table 2 

Overturning rates for AO ^ 1.7 K 
The same parameter values are assumed as in Table 1, 


p/p o w(cm s" 1 ) K(cm 2 s~ 1 ) t (s) T*/T 


0.1 

0.33 

6.6 

X 


6.1 

X 

10* 

3.2 

X 

10 

0.3 

0.11 

3.1 

X 

1°| 

2.6 

X 

10 L 

10 8 


X 

10 

1.0 

0.033 

1.3 

X 

10 5 

1.2 

X 

6.2 

X 

10' 


Notice that the time scale near the top is on the order of 
60 days, whereas in the mixing length convection case it is about 

Q 

an hour. Using a length scale L 2 x 10 cm, we estimate a 
meridional velocity of about 30 cm s”\ still smaller than the 
mixing length velocities even in spite of the large horizontal 
scale . 


, v 


* 

* 
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*0 The overturning time and the growth rate for the belts and 

o 

zones were estimated to be on the order of 10 s, while the 
horizontal diffusion time scale In deep layers was about half 
this. We therefore expect that the concentration of water vapor 
is Influenced by the belt-zone circulation to appreciable degree 
at levels just below the cloud base, but very little at levels 
many scale heights below. 


2k 


IV. Discussion 

The general picture we have outlined la that Jupiter's 
belts and zones represent convection driven by the Internal 
heat (and the fraction of the solar flux thnt penetrates to 
low levels). The phase change at the water cloud level pro- 
vides a peculiar kind of interface across which rising cur- 
rents undergo a temperature increment (relative to the dry 
adlabat). The thermal contrasts generated this way permit 
convection to c^.rry the necessary heat flux with very slow 
motions compared to those mixing length theory would predict. 
The planetary vortlcity gradient permits the shear associated 
with the axlsymmetric belts and zones to be stable; this ex- 
plains the symmetry and the length scale, as suggested by 
Ingersoll and Cuzzl ( 1969 ). 

Concerning motions of short time scales we have suggested 
that Inertial Instabilities probably occur and maintain the 
Richardson number larger than unity. We have argued that 
barotropic shear instability occurs, Bofcn these processes 
would lead to adjustments in the mean thermal field, and we 
showed that motions resulting from such adjustments would 
penetrate deep into the lower regions below the water cloud. 

It seems probable that when the sense of the vertical motion 
field is such as to draw moist gas upward, these sudden ad- 
justments could lead to an eruption of spots along an entire 
belt or zone, and these events are indeed observed. 

We wish now to ask what conditions are necessary for 
the existence of this kind of convection, and w!*ether they are 
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met on the other outer planets also. We begin by evaluating 
several numerical parameters. Assume that the potential 
temperature gradient la given .In the Intermediate layer by 

^ - k Ad sin ky , (34) 

where k Is the latitudinal wavenumber. Then 


—^•u(p 1 ) - ak 3 u Q sin ky , 


from (32) where 


u. 


Hae 


f o a 


(35) 


(36) 


is a thermal wind speed depending only on the planetary radius, 
a, and other physical parameters of the system. By substituting 
(35) into (33) we can solve for the critical wavenumber deter- 
mined by the barotropic instability criterion, and then for 
the maximum wind speed U. We obtain 

k3 * ^ , 0 - aku Q . (3?) 

We tabulate (37) for Jupiter, Saturn, and Uranus in 
Table 3 below. The corlolis parameters f and 3 are evaluated 
at 30 degrees latitude. The tropopause and water cloud pressures 
P x and p Q are based on the calculations presented by Weiden- 
schilling and Lewis (1973). The value of A0 is calculated from 

A6 * , (38) 

U P 
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where x* Is the water vapor number fraction, L is the latent 
heat (4.7 x HO 31 erg mole"" 1 ), m is the mean molecular weight. 


and Cp is the heat capacity at constant pressure. For a 


solar mixture, x’ * 1.0 x 10“^, m » 
erg K“ 1 gm” 1 , and we find A6 * 1.7K. 


2.3 Em mole 


-1 


1.3 x 10 


8 


We shall now calculate the ratio of vertical to horizontal 
potential temperature contrast, A0 v /A0, in the intermediate 
layer under the assumption that the Richardson number is unity, 
as suggested in Section II. The Richardson number is given by 


Ri 


fife* + JL\ 

T(3z + cl 


3U 

3z 


f 2 1 

- - R 

P 

K 

p 


a yj 


(39) 


where z measures geometrical height and the hydrostatic relation 
has been used to transform to the right hand expression. As- 
suming that 3 ®/3y is Independent of height (which is consistent 
only if A0 y << A0), and assuming that the maximum value of 3 0 /3y 
is kA0',* we can solve (39) for 3 ©/3p and integrate to obtain the 
maximum vertical potential temperature contrast. We obtain 


A0 


V _ k 2 RA0 1 -(t>i/ p 0 ) K 


A0 


** e 


(40) 


where e = kU/f. This verifies the estimate made in Section II, 
at least for one particular structure assumption. 

We see that as we move from Jupiter to Saturn to Uranus, 
assuming solar composition, both the wavelength-radius ratio 
and the Rossby number become larger. If these numbers approach 
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> 

1 

Table 3 

Dynamical parameters and cloud levels 


The first t';ree rows present data, based on Weidenschilling 
and Lewis (1973). H is the scale height at the water cloud 
temperature & . The next three rows present derived quantities 
as discussed in the text. The wavelength X * 2ir/k. Quantities 
in parenthesis for Saturn and Uranus are based on water abun- 
dances 2 and 5 times the solar composition value respectively. 


In all ( 

3ases, f 0 » 

-4 -1 

1.7 x 10 s and c » 1. 

r 

3 x 10 

8 erg K" 1 

gm” 1 . 


a(10 8 cm) 

—2 

g(cm s ) © 

H(km) p 

Q (atm) 

p l /p o 

g(cm”^s”^) 

Jupiter 

70 

2400 275 

43 

7 

14 

4.2 x lO" 1 ^ 

Saturn 

60 

1000 290 

107 

20 

40 

4.9 x 10" iq 

Uranus 

22 

900 360 

143 

150 

500 

1.3 x 10“ 13 


u o (ms" 1 ) 

U(ms -1 ) 

X(10 8 cm) 


X/a 

e 

Jupiter 

1.0 

28 

16 

.23 


.06 

Saturn 

1.4 (2.8) 

33 (52) 

16 (20) 

.27 

(.33) 

.09 (.10) 

Uranus 

4.5 (22) 

50 (150) 

12 (21) 

.56 

(.94) 

.15 (.26) 


unity we expect qualitatively different behavior for either 
geometrical or dynamical reasons. The principal parameter 
causing the increase is the decreasing planetary radius. In 
particular, X/a exceeds one-half for Uranus. (V/e have assumed, 
of course, that there are Internal heat sources or sufficiently 
deep penetration of sunlight to drive motions.) 
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Motions predicted on Saturn are only slightly stronger 
than on Jupiter, for solar composition. Observations suggest, 
however, that motions may b> several times stronger (Moore, 1939). 
The only explanation consistent with the Ideas proposed here Is 
that the 0:H ratio Is higher on Saturn. We have shown approximate 
values of derived parameters for Saturn and Uranus atmospheres 
enriched 2 and 5 times, respectively, relative to solar composition. 

Notice that because the latitudinal scale is determined 
by the barotropic Instability criterion, the last two columns 
in Table 3 are related. In fact, since 8 ^ ft/a, we have 
e 'v (ka)" 1 . 

Two comments can be made regarding the great red spot 
and the equatorial Jet. Our ideas here are consistent with 
Ingersoll’s (1973) suggestion that the red spot is a free 
atmospheric vortex. As he points out, once it is known 
that the barotropic vortlclty equation governs the flow, the 
question becomes one of initial conditions. We add one more 
question here: Can a vertical circulation pattern exist 

inside the spot to maintain the horizontal A9 necessary to 
provide the anti-cyclonic motion? 

The equatorial jet could be maintained by equatorward 
momentum transport in two-dimensional eddies governed by the 
same equation. Gierasch (1975) has made a similar suggestion 
regarding momentum transports in the Venus atmosphere. This 
question should be susceptible to numerical study. 

A number of other theoretical questions are raised by 
the conjectures 'we have made in this paper. One hesitates to 


make a long list when the observational basis is so meager. 

Two, however, are worth mentioning. We have assumed here that 
deep layers are horizontally uniform in composition and tempera- 
ture. What really .la the nnture of the coupling between the 
deep and intermediate layers? How do moisture and molecular 
weight affect small scale convection near the interface? Does 
a large scale circulation Indeed suppress moist convection? 

The second question is an easier one: Can the horizontal 

temperature gradients we postulate in the intermediate layer 
be stable? Does the presence of the deep layer suppress 
baroclinic Instabilities? 

Observational questions related to these ideas are 
obvious and we shall mention only a few. More details of the 
motion field at cloud level would be extremely useful, to 
determine whether small scales of motion (possibly related 
to inertial instabilities) exist, and to what accuracy the 
flow really is governed by vortlclty conservation. The thermal 
structure near the cloud tops is important. Composition 
variations (both water and other condensing constituents) may 
be extreme between zones and belts. Is the concentration 
greatest in zones? How deep do variations extend? Thermal 
contrasts in intermediate and deep layers are highly desirable 
but probably hopeless to obtain, since variations from adia- 
batic (both horizontally and vertically) are on the order 
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Figure Captions 


Figure 1, The three layers. An approximate temperature 
profile la indicated, with the moist adlabat 
sketched In as a dashed line. 

Figure 2. Development of belts and zones. In the Initial 
configuration convection occurs at all levels as 
indicated by the stippling. During growth, convection 
ceases near the water cloud level at locations of 
large scale downwelling. The final configuration 
has no small scale convection above the water cloud 
because the intermediate layer is stabilized as 
discussed in Section II. 
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